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An incident fast ion in the electronic stopping regime produces a track of excitations which can 
lead to particle ejection and cratering. Molecular Dynamics simulations of the evolution of the 
deposited energy were used to study the resulting crater morphology as a function of the excitation 
density in a cylindrical track for large angle of incidence with respect to the surface normal. Sur- 
prisingly, the overall behavior is shown to be similar to that seen in the experimental data for crater 
formation in polymers. However, the simulations give greater insight into the cratering process. The 
threshold for crater formation occurs when the excitation density approaches the cohesive energy 
density, and a crater rim is formed at about six times that energy density. The crater length scales 
roughly as the square root of the electronic stopping power, and the crater width and depth seem 
to saturate for the largest energy densities considered here. The number of ejected particles, the 
sputtering yield, is shown to be much smaller than simple estimates based on crater size unless the 
full crater morphology is considered. Therefore, crater size can not easily be used to estimate the 
sputtering yield. 

PACS numbers: 61.80.Jh, 79.20.Ap, 79.20.-m, 34.50.Fa 
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I. INTRODUCTION 



Surface modification of materials by 



irradiation has._b.een studied— ..in insulato 




and sur- 



semiconductors^-- and metalsE-Hi- using electron, STM 
and AFM microscopes. A large varie 
have been bservedu-'^«n|ms^tL craters 
crater rims (hillocks )Era_'Q_ii3 , _?l, adatoms__3, _ 
face roughening—^. Cratering occurs in response to the 
pressure pulse and fluid flow to the surface produced by 
the rapid deposition of energy, but the process is not un- 
derstood quantitatively. 'Bumps' generally appear when 
an energetic process occurs a few layers below the sur- 
face creating a low density region with a larger volume 
which raises the surface S3. When the energy loss per unit 
path length of the projectile, dE/dx, and the sputtering 
yield are relatively small, adatoms are observed in both 
experiments and simulations. By increasing the energy 
deposition (and the yield) craters are eventually formed. 
For very large energy deposition and yields, re-deposition 
of the ejecta plus plastic deformation occurs, producing 
craters with rims, studied recently for ion bombardment 
of polymers and other organic materials—. 

Craters are also produced by cluster ion bombardment 
which can lead to huge sputtering yieldsoEll. This pro- 
cess has been studied in the velocity regime in which 
nuclear (elastic) energy loss dominates over electronic en- 
ergy loss and-has been seen in both experimentsc3'E2l and 
simulationsOO'Ocj. The simulations are generally per- 
formed for bombardment at normal incidence and when 
energy is deposited in momentum transfer collisions to 



the target atoms. There ace-Jew simulations of crater- 
ing in the electronic regime— 3 and none for non-normal 
incidence. 

At normal incidence the crater produced by a fast in- 
cident ion has a roughly circular profile, but recent ex- 
periments have focused onions incident at a large angle 
with respect to the norma— I and at grazing incidence^. 
Even at normal incidence MD results for keV ion clus- 
ters incident on a copper surface appear to disagree 
with the scaling laws followed by macroscopic crateringllil 
in whieii the crater radius varies with the bombarding 
energyt- E as E\l 3 . For oblique incidence in polymers— 
and biomoleculesu bombarded by fast heavy ions it was 
found that the crater width does not increase significantly 
with increasing deposited energy density whereas the size 
of the crater along the incident ion direction increases 
rapidly with increasing energy density. 

In this paper molecular dynamics (MD) simulations 
are used to study the surface morphology produced by 
the energy deposited by fast ions incident at large an- 
gle with respect to the surface normal. The results of 
these simulations are compared to models for the length, 
width and depth of the crater vs. the energy density (i.e., 
dE/dx and track width). Since crater formation is used 
for sculptiag specific surfaces features for biopolecule 
adsoptionc _, for determining-iiiirface properties-, and for 
estimating sputtering yields— '0, we use MD simulations 
to extract scaling laws for crater formation. Although 
the simulations are for an "atomic" solid, quite remark- 
ably, the trends are very similar to those recently seen 
in polymers— However, in the simulations we can study 
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how cratering depends on the material properties and on 
the energy density deposited by an incident ion. 



II. MD SIMULATION 

Following the passage of a fast heavy ion a cylindri- 
cally energized region is produced in a solid, which we 
refer to as a track of excitations. A Lennard- Jones (L- 
J) crystalline solid is simulated with particles interacting 



through the potentials V (r) 

Although this is an oversimplified model of a real solid, 
this two parameter potential has the advantage that the 
equations of motion, and, hence, all results, including the 
crater dimensions, scale with e and a. In addition, cer- 
tain weakly bound solids, such as the low-temperature, 
condensed-gas solids, can be reasonably approximated as 
L-J solids, with parameters e and a taken to reproduce 
the material properties. All L-J samples have a cohesive 
energy U « 8e. The interlayer distance for (001) layers is 
d «s 0.78ct and the bulk modulus is B = 75e/a 3 
tails on the MD simulation can be found elsewhere^ 
As in our earlier papers, the scaling with s and a is re- 
placed by scaling using U and the number density n. For 
the fast processes which determine sputtering and cra- 
tering we showed that the scaling was roughly preserved 
when a more complex potential was usedE3. 

Since the results will be compared to data on polymers 
we note that certain polymers are roughly simulated us- 
ing a L-J potential for the inter-chain interactions, plus 
a stronger potential to account for the covalent interac- 
tion within-the chairO. A typical size for a in polymers 
is 3.5-5 AEj. The cohesive energy of a polymer is more 
difficult to define. The covalent bonds among atoms in 
the same polymer chain are of the order of several eV's 
with slightly weaker bonds between monomers in a chain. 
However, the bonding among atoms in the neighboring 
chains is very weak, much smaller than 1 eV, making 
the binding field "anisotropic" . Removing a small chain 
requires different energies depending on the chain orien- 
tation and entanglement. The average cohesive energy 
is usually taken tp-.be equal to the sublimation energy. 
A simple estimated gives U « 0.5 eV/monomer. There- 
fore, even though the L-J calculations scale with size and 
binding energy, we assume an effective binding energy of 
0.5 eV/particle and a =5 A so that our "atoms" very 
crudely represent monomers. The mass of the simulated 
particle only changes the time scale, which is given by 
the dimensionless time t/t , where t a = a^jM/s, and M 
is the mass of the simulated particle. Assuming a mass 
of 70u gives t a = 1.75 ps. 

The stopping power, i.e. the energy deposited per unit 
length dE/dx, and the track radius, here r cy i, are typ- 
ically used to describe the energy density deposited by 
the ion in its passage through the solid. Since in the 
electronic sputtering regime only a fraction of the exper- 
imental dE/dx goes into non-radiative de-excitations, in 



= 4e 



(r/cr) -(r/cr) 



the following we use the symbol (dE / dx) e ^ to represent 
the amount of energy deposition contributing to track 
formation, cratering and sputtering. This fraction has 
been estimated to be ~ 0.2 for MeV He + bombardment 
of solid O2E3. We also use this here to compare the exper- 
imental results with our MD simulations. To mimic the 
non-radiative energy release at the ion track in the MD 
simulations all N exc particles within a cylinder of radius 
r cy i are given an energy E exc with their velocities in ran- 
dom directions. Therefore, {dE/dx) e ^ is N exc E exc /d. 
A trac-k-,radius can be estimated from the Bohr adie 
radiusEZl, though this has been questioned recently^! 
Thermal spike models have been applied recently to esti- 
mate the latent-track radius in irradiated polymers with 
positive resultsa. In this paper all simulations were run 
for an initial track radius r cy i = la = 10 A, which implies 
Nexc ~ 10/ cos0. The incident angle, 0, is measured in 
degrees with respect to the surface normal. For all sim- 
ulations = 60. Larger angles were not feasible as the 
sample size required became too large to practically sim- 
ulate. Therefore, the stopping power from MD is multi- 
plied by a factor cos 60/ cos 79 to compare to experiments 
done at = 79. The size of the simulated sample was 
varied depending on the size of E exc such that the results 
did not depend on the boundary conditions used, and the 
total simulation time was also varied to be able to "de- 
tect" all ejected particles. Most quantities presented are 
averages of results from a number of simulations in which 
the directions of the energized atoms were randomly var- 
ied. 



III. CRATER FEATURES 



A cut-across a crater formed following an excitation 
event is shown in Fig. [l]. This cut is in a plane con- 
taining the initial surface normal and the track direction 
and shows the maximum depth of the crater. The crater 
wall is seen to have a slope similar to the incident ion 
direction on the entrance side and a very steep slope at 
the back. Remarkably, this shape is very stable even in 
this model solid as we have increased the run time by a 
factor of two and see little relaxation occurring. This is 
the case because of— tke rapid cooling in the track core 
discussed elsewhereSE]. At the end of a typical simu- 
lation (~ 30i o ) the temperature is ~50 K higher than 
the initial temperature of the sample (30 K). Of course, 
over very long times, relaxation can occur in an atomic 
material even at relatively low temperatures, whereas it 
is less likely in a polymer. Therefore, crater morphology 
in polymers can be readily studied. Despite some prob- 
lems with theLtinterpretation, scanning force microscopy 
measurementaij are regularly used to analyze heavy ion 
damage in solids. 
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FIG. 1. Slice of the final configuration of one MD cra- 
tering event, showing a side view of the crater and crater 
rim. The incident ion impact angle is O = 60, and 
(dE/dx) efl » 330U/a = 330 eV/nm. The different dimen- 
sions of the defect are shown in the figure. 

Since all craters were found to have a characteris- 
tic morphology, the dimensions used in the subsequent 
graphs and discussions are indicated in Fig. |l|. This 
crater is formed at high excitation densities (12nU) and 
exhibits a rim on the back side. Not all craters, however, 
have rims, as we will discuss. To 'see' thje rim better, in 
Fig. we show a top view of a cratercd, together with 
an AFM image of a heavy ion impacting on a PMMA 
film. The rim is primarily on the sides and the back of 
the crater, not at the ion entrance site. This is the case 
even when there is no momentum preferentially deposited 
along the incident angle as is the case in macroscopic cra- 
tering. Here the crater is formed from a cylindrical 'heat 
spike' but has many of the characteristics associated with 
impact cratering. We also note that atoms on the rim 
borders are aligned along the preferential [110] directions 
indicating recrystallization of the material pushed or de- 
posited onto the surface. For this "atomic" material a 
few adatoms are also seen far from the track region. 

The crater formation has several stages, but most of 
the crater volume is ejected before ~ 2Qt (35 ps for 
a mass of 70m). Temperature varies greatly during the 
formation process and near the center of the track it 
can be larger than the melting temperature even after 
10t o . The dependence of the crater dimensions on the 
energy deposited in the track can be seen in Fig. [|. 
The MD values represent the mean value of the crater 
dimension for 4-8 simulations at each (dE/dx) e jp and 
the maximum difference among the mean value and val- 
ues for particular simulation was taken as the error bar. 
Because of sample size limitations, angles above 60 de- 
grees are problematic whereas experimental results are 
often-performed at nearly grazing incidence. In previous 
workcj, the width of the distribution of original position 
of the ejecta along the direction of the incident beam was 
found to have a cos -1 dependence while no variation 
was found in the distribution along the perpendicular di- 
rection. Therefore, in order to compare the simulations 
with the experiments a cos -1 dependence has been as- 



sumed for lengths along the direction of incidence of the 
ion. As mentioned before, the MD results for crater and 
rim length obtained at = 60 were multiplied by a factor 
cos,60 / cos 79 when comparing to the data of Papaleo et 
alfl. Experimental results are the convolution of the ac- 
tual crater profile with the AFM tip shape, and may also 
involve some late relaxation of the crater walls. There- 
fore the measured depth is expected to be smaller than 
the actual depth. The effect of the convolution is smaller 
for the other crater dimensions. Finally, the initial track 
size is not only difficult to estimate but may vary with ion 
velocity. We have used r cy \ — 10 A in our simulations, 
but it is not unreasonable that the initial spike radius is 
5 times larger. For all these reasons the MD results were 
normalized separately for the length, width and depth 
comparisons in Fig 3. However, it is seen, quite remark- 
ably, that the trends in the experiment and in the 
simulation are the same. This indicates that useful 
scaling laws can be obtained and that the crater forma- 
tion process is insensitive to the details of both the energy 
deposition profile and the materials properties. This also 
means that simulations for relatively simple systems can 
be used to predict crater structures when using ions to 
modify materials. 



(a) 
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FIG. 2. Top view of a crater, where the color scale in- 
dicates height, (a) experimental result for O = 79, and 
(dE/dx) eff = 660 eV/nm (20 MeV Au on PMMA). (b) MD 
simulation for = 75 and (dE/dx) e ^ — 205 eV/nm. (c) 
MD simulation for 6 = 60 and (dE/dx).. = 330 eV/nm. 
Because of differences in B and energy, the MD crater in (c) 
is not as elongated as the experimental crater. 
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FIG. 3. Crater size as a function of {dE / dx) MD 
results for = 60, U = 0.5 eV, a = 5 A, and taking 
(dE/dx) e jj = 0.2 (dE/dx). MD crater length is multiplied 
by a factor (cos 60/ cos 79) to account for the different inci- 
dent angle in the simulation. Open symbols are experimental 
data from Papaleo et au. 
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FIG. 4. Same as Fig. 3, but showing rim size as a func- 
tion of (dE/dx) e j:j:. MD sizes are multiplied by 5 in order to 
compare trends with the experiment. 



For the Yamamura et alHaH shock model given in the 
appendix, the volume of ejecta is determined by the en- 
ergy density deposited, and the crater length and width 

are both expected to increase as ^J{dE/dx) e ^. Assum- 
ing that the volume removed and the crater size are di- 
rectly related, this gives a {dE/dxfJ^ dependence for 
the sputtering yield. The crater length is seen in Fig. 3 

to increase roughly as W 1 {dE/dx) e ^ in both the simula- 
tions and the experiments. However, the length appears 
to increase more slowly and to saturate at large values of 



(dE/da 



'eff 



After a steep increase, the crater width is 



seen to be much smaller than the length and increases 
only very slowly with increasing {dE/dx) e jf. Crater 
depth is much smaller than both the length and the width 
but appears to increase with {dE / 'dx) The pressure 
pulse model, also discussed in the appendix, gives such 
a scaling, but it also predicts that all of the dimensions 
have the same scaling. That model gave a good fit to MD 
calculations of the sputtering yield for ejection of large 
L-J molecules with hard cores, at normal incidence, and 
appeared to ragree with data from a solid made of large 
biomoleculesE3. However, it differs from what is seen in 
the simulations presented here and in the polymer exper- 
iments. 

Although it is not clear from Fig. 3, for E exc below 
U (i.e., the energy density in the track is less than the 
cohesive energy density) no crater is formed. Such a 
threshold is also seen in the experiments and, therefore, 
crater detection can give a measure of the cohesive en- 
ergy. Below the threshold, several atoms escape from the 
top layers leaving vacancies and the. sputtering yield is 
small, as discussed by Bringa et alSM. In addition, some 
atoms are displaced to the top layer, where they stay 
as adatoms, but an identifiable crater is not observed. 
For 6f7 > E exc > 2U a shallow crater forms and again 
several atoms are re-located as adatoms on the surface. 
For E exc > 6/7 the energy density is about a tenth of 
the bulk modulus of the material, slip dislocations ap- 
pear and a crater rim is formed. The dimensions of the 
rim are shown in Fig. [|as a function of {dE/dx) e ^. Af- 
ter the initial rise at 'threshold' the rim length and width 
stay constant within our error bars, but the rim height in- 
creases very slowly with {dE/dx) e jj. Again, these trends 
are also observed experimentally, which is quite remark- 
able considering the differences in materials. 



IV. SPUTTERING YIELD 



The sputtering yield can, of course, be obtained di- 
rectly from the MD simulations. We showed earlier the 
surprising result that at the high excitation densities for 
which the energized track produces craters, the sputter- 
ing yield, Ymd-, is not predicted by standard models. 
Ymd is roughly proportional to (dE/dx) e ^ times the 
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effective 'spjutter depth', which is a fraction of the initial 
track widtbO'lEl. When it is difficult to measure it di- 
rectly, the sputtering yield is often approximated by an 
estimate of the crater volume. This has been tested for 
normal incidence for L-J molecules with a corecll. Here 
we evaluate that procedure for large incident angle for 
a standard L-J solid. Typically one assumes that the 
ejected volume is a semi-ellipsoid, Y el1 = (ir/6) nri c r wc z c , 
where r c ;, r cw and z c are the crater length, width, and 
depth respectively. This rough estimate is based on 
crater shapes [Obtained from MD and has been used in, 
several papersB. In Fig. [| we show Y el1 for the polymersQ 
(Y£j. l p ) and the MD crater (Ymd) using the MD values for 
the crater in Fig. 3. These are both compared to the ac- 
tual MD yield, Ymd, as a function of (dE / dx) e j * . 



els for sputtering. However, it is seen that these yield 
estimates are an order of magnitude larger than Ymd, 
and that they have a steeper dependence on (dE/dx) e ^. 
This large discrepancy is due to several factors. First, 
the craters are not well approximated by a half ellipsoid. 
Second, for the excitation energies shown many of the 
atoms originally in the crater relocate on the rim. Third, 
there are regions of higher density at the crater walls. 
These factors add up to a surprisingly large overestimate 
of the experimental yield. More importantly, since the 
ratio Y^f D /Y MD changes with (dE/dx) e ^, the depen- 
dence of the yield with (dE/dx) e ^ can not be obtained 
from such estimates. Since most of the atoms on the rim 
come from the crater region, a more detailed description 
of the morphology of the crater needs to be made to ob- 
tain a reasonable yield estimate. 




(dE/dx) 8f) (eV/nm) 



FIG. 5. Sputtering yield as a function of (dE / dx) e ^. 
Yields obtained assuming the crater is an ellipsoid: Y^ v from 
experiments of ion bombardment of polymersu (solid squares) , 
Ymd fr° m the size of MD craters in Fig. 3 (open squares). 
The "true" MD yields, Ymd, are also shown (open circles). 
Lines are only a guide to the eye. 

The polymer yield was obtained as in Papaleo et aM 
using a semi-ellipsoidal volume and the same number 
density as in the MD simulations (using M = 70u gives a 
density of 0.98 g/cm 3 ). The experimental yield estimate 
for polymers is larger than the 'yield' estimated from the 
MD craters. The MD yields have been multiplied by a 
constant C, C — 2.5 x cos (60) / cos (79), which results 
from the normalization of crater size in Fig. 3. However, 
the trends are the same. They both give a yield that de- 
pends quadratically on (dE/dx) e ^. This is, fortuitously, 
the same dependence predicted by thermal spike mod- 



V. DISCUSSION AND SUMMARY 

Here we carried out a series of MD simulations to study 
crater formation due to the high energy density deposited 
in a cylindrical 'track'. Such an energized track might be 
formed by a penetrating fast ion that deposits its en- 
ergy in electronic excitations, which is of interest here, 
or deposits its energy by momentum transfer producing 



escribed are not im- 
, rather they are the 



recoil atoms. That is, the crater 
pact craters like the lunar craten 
craters formed in response to the rapid energy deposi- 
tion in the track of an energetic ion. The study here is 
for large angles of the track with respect to the surface 
normal. 

First, we showed that the crater structure remains sta- 
ble in this Lennard-Jones material over the longest sim- 
ulations time tested (~ 75i ), which is much larger than 
the crater formation time (~ 20t o ). This is the case be- 
cause of the rapid cooling of the track by the melting-^ 



pressure pulse processes described in earlier papersEJ 
We also found the initially surprising result that the scal- 
ing of the crater parameters with (dE / dx) e j * in this L-J 
solid agrees remarkably well with that found experimen- 
tally, for MeV heavy ion -bombardment of polymers at 
79 degrees to the normalfa This means that concepts 
learned from MD simulations of simple materials can be 
applied to more complex materials. 

As shown earlier for the experimental data for 
polymersB, we find here a threshold for crater formation 
and a second threshold for rim formation. In another 
set of experiments the rims could be removed when the 
polymer is maintained at higher temperatures so that vis- 
cous relaxation occurs. Therefore, Papaleo et aljj used 
the relaxation of rim formation vs. material tempera- 
ture to locate the glass transition temperature. Here we 
did not vary the material temperature as late relaxation 
occurring over long time periods can not be described 
using MD. However, in this paper we are able to relate 
the two thresholds to the track energy density. Our MD 
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simulations show that the threshold for crater formation 
occurs when the energy density in the track is close to 
the cohesive energy density or, in the track formation 
model used here, when the non-radiative relaxation en- 
ergy per particle inside the initial track, E exc , is near 
the sublimation energy U. The threshold for rim forma- 
tion, however, occurs at a higher energy density both in 
experiment and in the simulations. We find this to be, 
E exc ~ 6f, which occurs when the energy density in the 
track roughly equals the bulk modulus of the material. 
In the crystalline material slip dislocations can form at 
such energy densities allowing the raised structure to be 
maintained. Using an efficiency of 0.2 for determining 
(dE/dx) e ff and our simulation parameters, we find that 
the experimental value of the stopping power needed to 
form a rim would be around 2 keV/nm for = 79. 

Above the threshold the crater width is found to be 
nearly constant for large incident angles and the crater 
length and depth increase sub- linearly with (dE/dx) e j^. 
The rim height is ten times smaller than its length and 
grows faster than the crater size (depth and length). 
This dependence should be compared to the steeper 
dependence for the crater radius and depth at normal 
incidence£3. 

We note that for non-penetrating cluster bombard- 
ment the crater scaling with energy deposited is different 
from that found here. For incident clusters the energy 
of the projectile is deposited close to the surface. For a 
projectile energy E, the crater volume, V, is found to 
follower V oc (E/U). However, recent MD simulations 
seem to indicate that for keV copper clusters on copper 
scaling isE-3 V oc (E/U 2 ). It was argued that the pres- 
ence of a molten region in the MD simulation caused the 
steeper dependence on U . New MD results from simula- 
tions of keV xenon ions on golco support the quadratic 
dependence with U and relate this to the formation of a 
melt. It is difficult to compare these results with ours in 
which the energy is deposited in a long cylindrical track. 
Assuming that the energy relevant for crater formation 
is deposited in a volume close to the surface of depth L, 



prisingly large errors in the sputtering yield and, even, 
the wrong dependence on (dE/dx)*^. Therefore, if the 
yield can not be measured directly, the full morphology 
to the track and rim need to be described to get an ac- 
curate yield. 
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APPENDIX A: SHOCK MODELS FOR 
SPUTTERING AND CRATER FORMATION 



There are several related models which attempt to ex- 
plain sputtering at high excitation density using shock 
waves. Also, the collision of an impactor with a tar- 
get procbicing spallation has been extensively studied 
with MDc3, .together with cluster bombardment-induced 
shock wavesEj'tiZl. The spallation process originates from 
the interaction of two rarefaction waves, one coming from 
the shock wave reflectaeLat the surface and the other cojb-i 
ing from the impactorCJ. Yamamura and coworkersc3H 
estimated the sputtering yield due to shock waves with 
spherical symmetry intersecting a surface. They sug- 
gested that a hemispherical volume is ejected with 
radius r c . Then the yield, Y is proportional to the vol- 
ume of the ejecta (~ 27rrj?/3) with r c ~ (dE/dx)]/^ . 
Bitensky and ParilisE! considered cylindrical tracks and 
the incident angle dependence to model biomolecule sput- 
tering. At normal incidence their model reduced to the 
spherical shock model. The crater dimensions are as- 



E = (dE/dx) eff L, we find V oc (E/U) 2 . Here we also sumed to be proportional to yj (dE/dx) eff . 



find that the molten region is important in the crater for- 
mation. Clustering of the ejecta might also affect crater 
sizeEHl, but in our simulations, unlike the EAM Cu used 
in several cluster bombardment simulationaij, there is al- 
most no contribution of clusters to the sputtering yield. 

Finally, we examined the accuracy of roughly esti- 
mating the sputtering yield by simply parametrizing the 
crater volume. Recent results by Insepov et aZ.E3 point to 
a possible connection between crater size and hardness, 
and claim that crater volume is also related to the sput- 
tering yield. However, they find a different dependence 
on the bombarding energy for the yield (Y oc E 1A ) and 
for the crater size (V oc E), confirming the discrepancy 
found here (Y oc (dE/dx)^^ and V oc (dE/dx) e ^). We 
showed that, for the model material studied here, using 
crater size to estimate sputtering yield can produce sur- 



In order ta_explain experiments on ejection of whole 
biomoleculestZI where Y oc (dE/dx) e ff, Johnson et al£li 
proposed the pressure pulse model (PP). In the PP 
model, there are many excitation events along the ion 
track, each contributing to (dE/dx) e jf as in the simula- 
tions described here. Whereas the energy density evolves 
diffusively the net energy density gradient causes a net 
volume force and, therefore, a net momentum transfer 
radially and towards the surface. If the net momentum 
transfer to a certain volume is larger than some criti- 
cal momentum that volume will be ejected. This deter- 
mines a critical radius r c oc (dE/dx) e ff, with the vol- 
ume ejected proportional to r™. The PP model predicts 
an angular distribution peaked at 45° and agrees well 
with MD simulations that use a Lennard-Jones potential 
with a core to describe the interactions of large excited 
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molecules0. Notice that the PP model gives crater di- 
mensions proportional to AdE j 'dx) e * ^ . If a critical energy 
for ejection is considered!^, instead of a critical momen- 
tum, the yield is Y oc (dE / dx) 3 J^ , as in Kitazoe et aZ.@. 

In all models discussed above the width and length 
of the crater have the same dependence on (dE / dx) e ^ . 
When oblique ion incidence is considered the yield in- 
creases as 1/ cos 8 because the length of the crater in- 
creases also asE3 1/ cosG. Therefore, the (dE / dx) e ^ de- 
pendence is the same as at normal incidence, which is 
not what is found in the simulations presented here. 
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